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Abstract 

Individuals of different types, may it be genetic, cultural, or else, with different levels 
of fitness often compete for reproduction and survival. A fitter type generally has higher 
chances of disseminating their copies to other individuals. The fixation probability of 
a single mutant type introduced in a population of wild-type individuals quantifies how 
likely the mutant type spreads. How much the excess fitness of the mutant type increases 
its fixation probability, namely, the selection pressure, is important in assessing the impact 
of the introduced mutant. Previous studies mostly based on undirected and unweighted 
contact networks of individuals showed that the selection pressure depends on the struc- 
ture of networks and the rule of reproduction. Real networks underlying ecological and 
social interactions are usually directed or weighted. Here we examine how the selection 
pressure is modulated by directionality of interactions under several update rules. Our 
conclusions are twofold. First, directionality discounts the selection pressure for different 
networks and update rules. Second, given a network, the update rules in which death 
events precede reproduction events significantly decrease the selection pressure than the 
other rules. 
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1 Introduction 



In evolutionary dynamics, different types of individuals compete for survival in a population. 
A type means genotype, social behavior, cultural trait, and so on, depending on the context. A 
type that is fitter than others generally bears more offsprings. A major quantity representing 
how successfully a type spreads in evolutionary dynamics is the fixation probability (Moran, 
1958; Ewens, 2004; Lieberman et al., 2005; Nowak, 2006). In a simple case in which there are 
only two types, the fixation probability of a type is the probability that a single individual 
of that mutant type introduced in a population of the other wild-type individuals eventually 
occupies the entire population. Requirements for considering the fixation probability are that 
the evolutionary dynamics are stochastic and that the two unanimity states, that is, the one 
of the introduced type and the other of the wild type, are the only two absorbing states. The 
fixation probability of a type depends on the fitness of the type, connectivity of individuals, 
and the update rule of evolutionary dynamics (Ewens, 2004; Lieberman et al, 2005; Antal et 
ai, 2006; Nowak, 2006; Sood et al, 2008). 

Evolutionary dynamics, both ecological and social, pretty often occur on complex contact 
networks of individuals (Newman, 2003; Watts, 2004; Keehng 2005; Proulx et a/., 2005; May, 
2006). Some networks as well as the rule of reproduction and other factors amplify the selection 
pressure in the sense that a fitter type has a larger fixation probability and a less fit type has a 
smaller fixation probability compared to a reference case of the all-to-all connected population 
(Lieberman et al, 2005; Antal et al, 2006; Sood et al, 2008). Other combinations of a network 
and an update rule may suppress evolutionary pressure, with the fixation probability relatively 
insensitive to the fitness of the mutant type. To quantify the extent to which a particular 
situation amplifies or suppresses the selection pressure is important for assessing the impact of 
a mutant type. 

Studies of fixation probability on networks and structured populations have been restricted 

to two neutral types (Donnelly and Welsh, 1983; Taylor, 1990; Taylor, 1996) or undirected and 
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unweighted networks (Mamyama, 1970; Slatkin, 1981; Antal et ai, 2006; Sood et ai, 2008), 
albeit a notable exception (Lieberman et ai, 2005). Regarding the directionality, undirected 
(and unweighted) networks may be natural for modeling social interaction such as games, where 
an adjacent pair of individuals, for example, is simultaneously involved in a single game to 
determine the fitness of each individual. However, reproduction events given the fitness of each 
individual may occur on directed or weighted networks. Indeed, relevant contacts in many real 
networks have directionality because of, for example, heterogeneity in the size of habitat patches 
(Gustafson and Gardner, 1996) and geographical biases such as the wind direction (Schooley 
and Wiens, 2003) and riverine streams (Schick and Lindley, 2007) in ecological networks. Social 
networks based on glooming behavior of rhesus monkeys (Sade, 1972), email communication of 
humans (Ebel et ai, 2002; Newman et ai, 2002) also have directionality. In the present work, 
we examine the evolutionary dynamics on different undirected and directed networks, namely, 
the complete graph, the undirected cycle, the directed cycle, the weighted undirected star, 
undirected and directed random graphs, and undirected and directed scale-free (SF) networks. 
Motivated by previous studies (Antal et ai, 2006; Ohtsuki et ai, 2006; Ohtsuki and Nowak, 
2006; Sood et ai, 2008), we examine the effects of several update rules on the selection pressure. 
We show that asymmetric connectivity generally turns down the selection pressure and that 
specific update rules suppress the selection pressure more than other rules. 

2 Model 

2.1 Dynamics 

We consider a population of N haploid, asexually reproducing individuals. The structure of 
the population is described by a directed graph G. We denote by E the set of directed edges, 
implying that {vi,Vj) & E ii and only if there is a directed edge from Vi to vj. Each node Vi 
(1 < i < N) is occupied by an individual of one type. Edges represent the likelihood with 
which the types are transferred from nodes to nodes in evolutionary dynamics. 
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We assume that there are two types A and B. A node is occupied by either type A or B. 
Types A (mutant type) and B (wild type) confer fitness r and 1 to its bearer, respectively. 
Given graph G and an update rule, the collective network state, which is specified by an 
assignment of A or i? to each node of the network, stochastically evolves. The dynamics last 
until the unanimity (i.e., fixation) of A or that of B is reached. We do not assume mutation 
so that these two collective states are the only absorbing states. We are concerned with the 
fixation probability of type A denoted by Fy{r). It is the probability that a single type- A 
mutant introduced at node f in a sea of — 1 resident type-i? individuals fixates. Type B 
fixates with probability 1 — -F^(r). The fixation probability depends on G, the update rule, and 
the initial location v of the A mutant. To focus on the effect of G and the update rule, we 
examine F{r) = Fy{r)/N . The effect of the initial location of the mutant has been analyzed 
for undirected (Antal et al, 2006; Sood et al, 2008) and directed (Masuda and Ohtsuki, 2009) 
networks. 

An update event occurs on one directed edge per unit time. We assume that the direction 
of reproduction is the same as the edge direction. We examine the following five update rules, 
most of which are motivated by past literature. For the sake of explanation, we explain the 
update rules for directed and unweighted networks. However, we can extend the model and 
the results to the case of weighted edges in a straightforward manner (Lieberman et al, 2005) 
(also see related analysis in Sec. 14. 4p . 

2.2 BD-B update rule 

Under the birth-death rule with selection on the birth (BD-B), we first select one node Vi for 
reproduction in each time step. The probability that Vi is selected is proportional to the fitness 
value, that is, fi/J2i fh where fi G {r, 1} is the fitness of the type on Vi. Note that selection 
operates on the birth. Then, the type at Vi is propagated to a neighbor of Vi along a directed 
edge that is chosen with probability 1/A;°"*, where /c°"* is the outdegree of fj. The probability 
that directed edge (fj, Vj) & E is used for reproduction is equal to fi/k°^^ Ylii fi- 



5 



BD-B is the update rule considered in (Lieberman et ai, 2005). It is equivalent to the 
previously defined birth-death rule for games on undirected networks (Ohtsuki et ai, 2006; 
Ohtsuki and Nowak, 2006) and to the invasion process defined in (Antal et al., 2006; Sood et 
ai, 2008). 

2.3 BD-D update rule 

Under the birth-death rule with selection on the death (BD-D), a random individual Vi is first 
chosen for reproduction with equal probability 1/N. Then, one of its neighbors Vj that receives 
a directed edge from Vi dies with probability proportional to 1/ fj, and the type at Vi replaces 
that at Vj. Selection operates on the death. The probability that edge {vi,Vj) e E is used for 
reproduction is equal to /r^- When r = 1, BD-B and BD-D are identical. 

2.4 DB-B update rule 

Under the death-birth rule with selection on the birth (DB-B), a random individual Vj first 
dies with equal probability 1/N. Then, a neighbor of vj that sends a directed edge to vj, 
denoted by Vi, is selected for reproduction with probability fi/ J2i-{vi vj)eE fi' ^^'^ ^JP^ of 
Vi replaces that at Vj. Selection operates on the birth. The probability that {vi,Vj) & E is 
used for reproduction is equal to fi/NY^^.^^^ vj)eE fi- DB-B is equivalent to the score-dependent 
fertility model proposed in (Nakamaru et ai, 1998) and to the death-birth rule previously 
used for evolutionary games on undirected networks (Ohtsuki et al, 2006; Ohtsuki and Nowak, 
2006). Many numerical studies of spatial reciprocity (Nowak and May, 1992) and network 
reciprocity (Santos and Pacheco, 2005) are also based on this update rule or similar rules. 

2.5 DB-D update rule 

Under the death-birth rule with selection on the death (DB-D), a node Vj is first chosen for 
death with probability fj^/^ifj^^- Selection operates on the death. Next, a neighbor of 
Vj that sends a directed edge to Vj, denoted by Vi, is chosen for reproduction randomly with 
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probability where kj^ is the indegree of Vj. Then, the type at Vi replaces that at Vj. The 

probability that {vi,Vj) E E is used for reproduction is equal to fj'^/kj^'^iff^- This rule is 
equivalent to the score-dependent viability model proposed in (Nakamaru et al, 1997) and to 
the voter model defined in (Antal et al, 2006; Sood et al, 2008). When r = 1, DB-B and 
DB-D are identical. 

2.6 LD update rule 

Under the link dynamics (LD) with selection on the birth (resp. death), directed edge (vj, Vj) e 
E is chosen with probability h/ Y.{vi,vi,)eE fi (^esp. //V E(^;,,„,)6B /r^)- The denominators in 
these probabilities indicate the summation over all the directed edges. Selection operates on 
the birth (resp. death). Then, the type at Vi replaces that at Vj. When there are only two 
types, as assumed in this and most previous studies (Lieberman et al, 2005; Antal et al, 2006; 
Sood et al, 2008), the LD with selection on the birth and that on the death coincide with each 
other, up to a change of the timescale. Then this rule is the same as the LD defined in (Antal 
et al, 2006; Sood et al, 2008). 

3 Evolutionary amplifiers and suppressors 

In this section, we explain the order parameter measured in the following analysis. We are 
concerned with dependence of selection pressure on networks and update rules. The fixation 
probability at neutrality, namely, r = 1, is equal to F{1) — 1/N. A type with a larger 
fitness is more hkely to fixate, that is, dF{r)/dr > 0. However, how F(r) depends on r 
differs by the network and by the update rule. The reference profile of F(r) is given by the 
Moran process (Moran, 1958; Ewens, 2004; Nowak, 2006). In the Moran process, we select an 
individual Vi at each time step for reproduction with the probability proportional to /j. Then, 
the offspring replaces an individual randomly picked from the rest of the population with the 
equal probabihty 1/(A^ — 1). The evolutionary dynamics on the complete graph under BD-B is 
the Moran process. As explained in Appendix B, F{r) for the Moran process is equal to F{r) 
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for the complete graph under DB-D and LD as well as under BD-B, and it is given by 

1 -r-i 



(1) 



1 

A network or an update rule that yields F{r) larger than Eq. ([T]) for r > 1 and smaller than 
Eq. ([T]) for r < 1 is called evolutionary amplifier. For an evolutionary amplifier, the selection 
pressure is magnified relative to that of the Moran process (Lieberman et al, 2005). If the 
fixation probability is smaller than Eq. ([T]) for r > 1 and larger than Eq. ([T]) for r < 1, the 
network or the update rule is called evolutionary suppressor. 

The amplification factor K may be defined when the fixation probability is given in the 
form: 

1 - r"^ 

(2) 



I _ J.-KN' 

If Eq. holds for a graph with nodes, the evolutionary dynamics are equivalent to the 
Moran process with the effective population size KN . Because a larger population size results in 
stronger selection pressure for the Moran process, K > 1 and K <1 correspond to evolutionary 
amplifier and suppressor, respectively. For example, for a family of star graphs under BD- 
B (Lieberman et al., 2005), and for undirected uncorrelated graphs under BD-B and DB- 
D (Antal et al., 2006; Sood et al., 2008), F{r) obeys Eq. ([2]) with different values of K in an 
appropriate limit. However, F{r) generally deviates from Eq. ([2]). Therefore, in the numerical 
simulations in Sec. [5l we calculate the fixation probability at r = 4 and compare it with that 
of the Moran process. 

4 Fixation probability for some simple graphs 

For some simple graphs, we calculate the fixation probability for different update rules to 
determine the selection pressure. The details of the calculations are shown in Appendix B. The 
results in this section are mostly restricted to undirected networks. We will treat the effect of 
directionality numerically in Sec. O 
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4.1 Complete graph 

Consider the evolutionary dynamics on the complete graph of nodes depicted in Fig. [U^a). 
Self loops are excluded. The complete graph is undirected and regular, where a regular graph 
is one in which the indegree and the outdegree of all the nodes are the same. For undirected 
regular graphs, BD-B, DB-D, and LD are all equivalent to the Moran process with the same 
population size N (Antal et al, 2006; Sood et al, 2008). Therefore, F(r) is given by Eq. ([1]). 
For BD-D, we obtain 

^^^^ = 1 I (AT ^'^^ 

± -t- I^JV L) Z^m=l rn+r{N-l-m) 

For DB-B, we obtain 

In Fig. [21(a), F{r) for the Moran process, BD-D, and DB-B is compared for the complete 
graph with = 10. The three lines cross at neutrality, that is, (r, F(r)) = (1,1/A^). The 
complete graph is an evolutionary suppressor under BD-D, even though lim^^oo F{r) = 1 
and lim^^o F{r) = 0. This is because selection occurs among the — 1 nodes excluding the 
reproducing node. DB-B is a stronger suppressor than BD-D. Differently from BD-D, we obtain 
lim^^oo F{r) = {N — 1)/N and limr_+o F{r) = for DB-B, so that a mutant may not fix under 
DB-B even if its fitness is infinitely large. 

For the three different cases, F{4) is plotted against A^(> 2) in Fig.[2](b). The fixation prob- 
ability for BD-D and DB-B converges to that for the Moran process as A^ ^ oo. Remarkably, 
for BD-D and DB-B, the existence of more competitors in the population, that is, larger A^, 
leads to the higher fixation probability of the single type- A mutant. As A^ increases, selection 
acts on more nodes relative to the population size (i.e., (A^ — 1)/A^) under BD-D and DB-B. 
Then tjpe-A nodes, which are rare in an early stage of dynamics, are involved in competition 
for survival (under BD-D) or reproduction (under DB-B) more often such that type A takes 
advantage of being inherently fitter than type B. 
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4.2 Undirected cycle 

Consider the undirected cycle of size N depicted in Fig. [U^b). Because the undirected cycle is 
unweighted and regular, BD-B, DB-D, and LD are again equivalent to the Moran process, and 
F{r) is given by Eq. ([T]). For BD-D, we obtain 

1 - r"^ 



For DB-B, we obtain 



1 , r-l , (r2-2r-l)r-^+l 
r(r+l) r+1 



1 — r" 



= 1 I r-l I (.-3)r-^ - - (6) 

2'r 2 

For the undirected cycle, F(r) with = 10 and -F(4) with varied are shown in Fig. [2t^c) 
and Fig. [2](d), respectively. Qualitatively agreeing with the case of the complete graph, BD-D 
is suppressing, and DB-B is even more so. In contrast to the case of the complete graph, BD-D 
and DB-B persist to be suppressing for large iV. Under BD-D and DB-B, selection operates on 
about (A;) individuals, where (fc) is the mean degree of the network, whereas selection operates 
on N individuals under BD-B, DB-D, and LD. For the complete graph, the difference diminishes 
as — ^ oo because (k) = N — 1. For the undirected cycle, (k) = 2 independent of A^, which is 
a likely reason why BD-D and DB-B are suppressing even for large A^. 

4.3 Directed cycle 

Consider the directed cycle of size A^ depicted in Fig. [T](c) . It is straightforward to verify that 
the evolutionary dynamics under BD-B, DB-D, and LD are equivalent to the Moran process. 
For BD-D and DB-B, selection pressure is totally annihilated, that is, F{r) = 1/N. 

4.4 Star 

Consider the star with A^ nodes depicted in Fig.lH^d). One central hub is connected to the other 
A^ — 1 leaves. Only in this section, we introduce the edge weight for a computation purpose. 
Specifically, each edge outgoing from the hub has weight 1, and each edge incoming to the hub 
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has weight a. The edge weight is assumed to be muhiphed to the hkehhood with which the 
edge is used for reproduction (see Sees. I2.2ti2.6l and Appendix B). 

The fixation probabihty for the weighted star under LD is given by 



N-l r+a 1 ra+l 



1 _ r-N ( 1±3l\ 
^ ' \ra+l) 

For general a, we obtain in the hmit N ^ oo 



In 



r+g 



K = (8) 

mr 

where K is defined by Eq. ([2]) . Although in an incomplete form of Eq. ([H]) , which includes r in 
the RHS, K > 1 {K < 1) for all r 7^ 1 when a > 1 (a < 1). Therefore, the weighted star is an 
amplifier (a suppressor) when a > 1 (a < 1). With a = 1 in Eq. ([7]), F(r) is equal to Eq. ([1]). 
This is expected because the evolutionary dynamics for any undirected network are equivalent 
to the Moran process under LD (Antal et ai, 2006; Sood et ai, 2008). 

As shown in Appendix B, LD on the weighted star with a = A^ — 1 is equivalent to BD-B on 
the unweighted star. Therefore, the unweighted star is an amplifier under BD-B. Particularly, 
in the limit A^ — 00, Eq. yields K = 2, consistent with the previous result (Lieberman et 
ai, 2005). The theory for undirected uncorrelated networks (Antal et ai, 2006; Sood et ai, 
2008) also predicts K = (k) {k~^) ^ 2 for the undirected star under BD-B. 

Also shown in Appendix B is that LD on the weighted star with a = 1/ {N — 1) is equivalent 
to DB-D on the unweighted star. Setting a = 1/{N — 1) in Eq. ([8]) yields = in the 
limit N 00. This is consistent with the result for undirected uncorrelated networks: K = 
{kf I ik^) = 2/N (Antal et ai, 2006; Sood et ai, 2008). 

The fixation probability for BD-D is given by 

Fir) = -— V, (9) 



N 1 I l+r(Af-l) (. _ r-N+^N 



and that for DB-B is given by 



_ {rN + N + 2r - 2){rN - r + 1) 

~ m{N + 2r-2){r + l) ' ^ ' 
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For the unweighted star, F{r) with = 10 and F{A) with varied are shown in Fig. M^e) 
and Fig. [21(f), respectively. LD corresponds to the Moran process, and BD-B is the only 
amplifying update rule. Among the other three suppressing update rules, DB-B is the most 
suppressing. Figure [2](f) indicates that BD-D is more suppressing than DB-B for small 
and vice versa for large A^. As A^ tends large, DB-B and DB-D become strongly suppressing. 
Indeed, both Eq. ([71) with a = 1/{N — 1) representing DB-D and Eq. ( |T0|) representing DB-B 
yield F{r) = 0{1/N), r > 1. In contrast, the fixation probability for BD-D approaches that 
for the Moran process as A^ — * oo; Eq. ([9]) yields limAr^oo -^(^) = {r — l)/r, which agrees with 
the result for the Moran process derived from Eq. ([T]). 

5 Numerical results 

Here we report numerical results for the fixation probability of the mutant type with r = 4. 
Equation ([T]) implies limTv^oo -^('") = 3/4 for the Moran process. If F{A) measured for a 
combination of a network and an update rule is larger (smaller) than 3/4, that combination 
is probably amplifying (suppressing). Whether the combination is amplifying or suppressing 
and to what extent actually depend on r. However, our extensive numerical results suggest 
that this dependence is not strong unless r is extremely small or large. The values of F{4) 
are also well correlated with dF{r)/dr\r=i, which is the sensitivity of the selection pressure at 
neutrality. Therefore, we assume r = 4 in the following. 

A graph G is called strongly connected if there is a directed path from an arbitrary chosen 
node to another. For G that is not strongly connected, the mutant type introduced at a node 
in a downstream component never fixates, and the fixation problem is ascribed to that for the 
most upstream strongly connected component of G (Lieberman et ai, 2005). Therefore, we 
assume that G is strongly connected without loss of generality. 
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5.1 Small networks 

In this section, we present numerical results for small networks. The evolutionary dynamics 
are interpreted as a discrete-time Markov chain on a finite space. A state of the Markov chain 
is specified by an assignment of type A or 5 to each node, so that there are 2^ possible states. 
Among them, the two states corresponding to the unanimity of A and that of B are the unique 
absorbing states. Given G and an update rule, we obtain the exact fixation probability by 
solving a system of linear equations of 2^ — 2 variables (see Appendix A for methods). Because 
of the computation time, we set = 6 and calculate F{A) for all the 1047008 strongly connected 
networks under the five update rules. The results are qualitatively similar for A^ = 4 and A^ = 5 
(data not shown). 

To visualize the results and to correlate -F(4) with the structure of G, we regard the values of 
F(4) for different networks as different data points. Then we regress F{A) of these data points 
against various measurements of the network, which we call order parameters. We look for the 
order parameters that are more correlated with -F(4) than other order parameters. Our choice 
of the order parameters is arbitrary. For each of the five update rules, we list in Tab. [U the 
Pearson correlation coefficient between -F(4) and the order parameters, where (■) denotes the 
average over the nodes. The first order parameter is the mean degree {k) = (A;*") = (A;""*), where 
k^^ and is the indegree and the outdegree of a node, respectively. Seven order parameters 
are normalized moments of the degree distribution: (A;) (1/fc™), (A;) (l/A;""*), (A;)V (('^*")^), 
(A;)^ / {^(A;°"*)^y (A;)^ / (A;*"A;°"*), (A;°"VA;^'^), and (A;*'^/A;°"*). The other two order parameters are 
the normalized standard deviation of the node temperature: std{T^^) and std{T°^^). The tem- 
perature T?" of node Vi (Lieberman et al, 2005) is defined by T?" = J2j-{v- Vi)eE i'^^ji/ 
We define strf(T^") = ^Eii - {T^-'')f / N / (T^") = ^Eti iV-'^f /N, where we have 
used (T*") = 1. For isothermal networks, which satisfy std{T'^^) = 0, the evolutionary dynam- 
ics under BD-B are equivalent to the Moran process (Lieberman et al, 2005). Similarly, we 
define the temperature for the outdegree by T°"* = J2j-{v,vj)eE i'^^^j/ which satisfies 

(T°"*) = 1, and measure std{T°'^^). 
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Under BD-B, we have not found an order parameter that is strongly correlated with -F(4), as 
listed in Tab. [TJ For a visualization purpose, we plot F{4) against std{Tin) for all the networks 
in Fig. [3](a). Each gray dot in the figure corresponds to a network of size = 6. Black 
circles are for undirected networks. The Moran process with = 6 yields -F(4) = 0.7502, 
which is shown by the solid line. Because r = 4 > 1, networks with F{4) above (below) the 
solid line are amplifiers (suppressors). Figure [3]^a) suggests the following. First, the undirected 
star (indicated by the arrow) is by far the most amplifying among all the networks. Most 
networks that yield large F{4) are variants of the star and have large std{Tin). Second, all 
the undirected networks are amplifiers under BD-B, consistent with the result that F{r) for 
undirected uncorrelated networks is given by Eq. ([2]) with K = (k) {k"^) > 1 (Antal et ai, 
2006; Sood et ai, 2008). Indeed, F{A) for undirected networks with = 6 (black circles) are 
strongly correlated with (k) {k"-^) (data not shown). Third, a majority of directed networks is 
suppressor. The mean and the standard deviation of F{A) based on all the networks and those 
based on the undirected networks compared in Tab. [2] indicate a significant difference. 

Under BD-D, (k) is most strongly correlated with F{4). The values of F{4) for different 
networks are plotted against (k) in Fig. [S](b). The complete graph, which corresponds to 
(A;) = A^ — 1 = 5, is the least suppressing, although it is nevertheless a suppressor, as shown in 
Sec. 14.11 All the networks are suppressors. Table [2] indicates that the suppression is generally 
stronger under BD-D than under BD-B and that directed networks are stronger suppressors 
than undirected networks on an average. 

Under DB-B also, (k) is most strongly correlated with -F(4). Figure [3](c), which shows the 
relation between F{A) and (k), indicates that the complete graph achieves the largest F{4), 
that all the networks are suppressors, and that directed networks are more suppressing than 
undirected networks on an average. These tendencies are similar to those for BD-D. 

Under DB-D, std{T°^^) is most strongly correlated with F{4). Equation ([2]) with K = 
{ky /{e), which was derived for undirected uncorrelated networks (Antal et ai, 2006; Sood et 
ai, 2008), approximates our numerical results for the undirected networks (data not shown). 
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However, std{T°^*) seems to be a better predictor of F{A) than (k)"^ / ({k'^^)j, {k)^ / (|(A;°"*)^y, 
and (A;)^ / probably because of the small network size. Figure [3]^d) indicates that -F(4) 

for the networks with std{T°^^) = is equal to F{A) for the Moran process, that all the other 
networks are suppressors, and that directed networks are more suppressing than undirected 
networks on average. 

Under LD, is most strongly correlated with -F(4). Evolutionary dynamics for 

undirected networks are equivalent to the Moran process (Antal et ai, 2006; Sood et ai, 2008). 
Figure O^e) indicates that directed networks are generally more suppressing than undirected 
networks, that there are both amplifiers and suppressors among directed networks with the 
majority being suppressors, and that smaller tends to yield larger F{4). 

The network that achieves the largest F{4) is shown in Fig. HI Most networks that pro- 
duce large -F(4) are variants of this network. Because Xlili ^i" ~ YliLi k"^^, a small value of 
{k^^/k°'^'^) requires that A;*" < holds for many small-degree nodes and that A;*" > A;°"* holds 
for a relatively small number of hubs. The network shown in Fig. H] complies with this property. 
In Fig. m each of the two hubs links to the half of the peripheral nodes, whereas each peripheral 
node links to both hubs. If we merge the two hubs into one as an approximation, the network 
is regarded as a weighted star with a = 2 (Fig. [Hd)). The weighted star with a = 2 is indeed 
an amplifier, as shown in Sec. 14. 4[ 

Based on the numerical results for the five update rules, we claim the following. First, 
directed networks tend to be suppressors compared to undirected networks regardless of the 
update rule. Second, some update rules (i.e., BD-D, DB-B, and DB-D) are much more sup- 
pressing than the others (i.e., BD-B and LD). Particularly, the magnitude of the amplification 
is ordered as: BD-B > LD > DB-D > BD-D > DB-B (Tab. which is consistent with the one 
for the undirected star with small N (Fig. [21(f))- In a single time step, the selection pressure 
operates on nodes for BD-B, DB-D, and LD. However, it operates on at most — 1 nodes 
for BD-D and DB-B, because the node first selected for reproduction under BD-D or for death 
under DB-B does not participate in the competition. This is a main reason why BD-D and 
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DB-B are strongly suppressing. We will corroborate these points by numerical simulations of 
large networks in Sec. 15. 2[ 

5.2 Large networks 

Because of the computational cost, the exact numerical analysis performed in Sec. 15. H is feasible 
only for small networks. In this section, we examine evolutionary dynamics on larger networks 
by Monte Carlo simulations. 

We use the Erdos-Renyi (ER) random graph, SF networks, and their directed versions of 
different sizes (e.g. Albert and Barabasi, 2002; Newman, 2003). We generate the undirected 
ER graph of mean degree {k) by connecting each pair of nodes with probability {k) /{N — 1). 
The directed ER graph is generated by connecting each ordered pair of nodes with probability 
2 (k) /{N — 1) so that (fc*") = = (k). In both cases, the degree distribution p{k) follows 

a poisson distribution with mean (k): p{k) ^ e"^/^'^^/ (k). For SF networks, we assume p{k;) oc 
k~^, {{k) /2 < k), where the power-law exponent 3 is an arbitrary choice. The SF network 
represents the situation in which the degree is strongly heterogeneous (Albert and Barabasi, 
2002; Newman, 2003). To generate a SF network, we first determine the degree of each node 
stochastically according to the power-law distribution p{k) with the restriction that the sum of 
the degree is even. Then, we randomly add edges one by one so that the predetermined degrees 
are respected at each node (Albert and Barabasi, 2002; Newman, 2003). For undirected and 
directed SF networks, the added edges are undirected and directed, respectively. In each of the 
four network models, we discard networks that are not strongly connected. 

To calculate -F(4), we perform 2000 runs for each initial location of the type A individual 
and count the fraction of the 2000iV runs in which the unanimity of type A is reached. This is 
a Monte Carlo realization of the fixation probability. For a fixed A^, (k), network model, and 
update rule, we calculate the average and the standard deviation of -F(4) based on 20 samples 
of networks. 

For (k) = 10, -F(4) under BD-B is plotted against A^ for the ER and SF networks in Figs. [!](&) 
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and (b), respectively. The standard deviation is shown by the error bars. The SF networks are 
more amphfying than the ER networks, consistent with the previous results (Lieberman et ai, 
2005; Antal et ai, 2006; Sood et ai, 2008). In addition, undirected networks (solid lines) are 
more amplifying than directed networks (dashed lines), even for large A^. The same tendencies 
are found for BD-D (Figs. El^c) and (d)), DB-B (Figs.[5](e) and (f)), DB-D (Figs. EKg) and (h)), 
and LD (Figs. [5](i) and (j)). This is a strong indication that directionality of edges generally 
suppresses the selection pressure. 

For each of the four types of networks used in Fig. [5l the strength of the selection pressure 
is ordered as: BD-B > BD-D, LD > DB-D > DB-B. This order is consistent with the one for 
the unweighted star with large (Fig. [2](f)). Compared to the results for the small networks 
shown in Sec. 15.11 BD-D is much less suppressing in large networks, as shown in Figs. [5](c) and 
(d). This is a bit surprising because we have kept (k) = 10 in Fig. [5] so that the competition 
for survival happens among roughly (k) <^ nodes per unit time. In contrast, DB-B remains 
strongly suppressing even for large A^ (Figs. [S](e) and (f)). An increase in (k) with A^ fixed 
makes DB-B less suppressing, as shown in Fig. [HI This is consistent with the behavior of the 
ensemble of the small networks (Fig. [31(c)). 

6 A theoretical explanation of the effect of directionality 

The finding that the directionality of the network suppresses the selection pressure can be 
explained for LD on a simple network. Consider a bidirectionally but asymmetrically connected 
two-node network. Edge (fi, V2) (from vi to V2) has weight unity, and edge (f2, fi) has weight 
a. Even though we have formulated the evolutionary dynamics on unweighted networks, the 
extension to the case of weighted networks is straightforward (also refer to the analysis of the 
weighted star in Appendix B). By enumerating the possible reproduction events per unit time. 
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we obtain 



F{AB) 



rF(AA)+aF(BB) 



r 



(11) 



F{BA) 



raF(AA)+F{BB) 
ra+1 



r+a 



r + a' 
ra 



(12) 



ra + 1 ' 



where the first and the second arguments of F are the types at Vi and V2, respectively, and 
F is the fixation probabihty starting from that network state. Note that F{AA) = 1 and 
F{BB) = 0. Using Eqs. ([II]) and (pSD, we obtain 



For r > 1, F{r) as a function of a takes the maximum at a = 1. For r < 1, F{r) takes 
the minimum at a = 1. Therefore, the unweighted network (i.e., a = 1) is the most amphfying 
when we vary a. LD on this unweighted network (i.e., a = 1) is equivalent to the Moran process 
with N = 2. Accordingly, weighted networks (i.e., a 7^ 1) are suppressors, and the suppression 
is stronger as a deviates from unity. 

The applicability of the arguments above is beyond the two-node network. Suppose that 
a network is divided into two modules with homogeneous intramodular connectivity and that 
intermodular connectivity is sparse and homogeneous, possibly with more edges from one mod- 
ule to the other than the converse. Then the fixation probability of a mutant on a particular 
node will depend only on the module in which the initial mutant invades. In this case, we can 
approximate the network by the two-node weighted network, where the edge weights between 
the two aggregated nodes represent the gross connectivity between the two modules. 



Motivated by the observation that many real ecological and social networks underlying repro- 
duction in evolutionary dynamics are directed or weighted, we have investigated the fixation 
probability on networks with directed edges under five update rules. For undirected networks, 
selection pressure relies on the network and the update rule, which reproduces the previous 




(13) 



7 Discussion 
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results (Lieberman et ai, 2005; Antal et ai, 2006; Sood et ai, 2008). Our main conclusions in 
the present paper are twofold. 

First, directionality of edges suppresses the selection pressure for all the network types and 
the update rules that we have dealt with. We have presented numerical evidence for different 
sizes of networks (Sec. [5]) and simple analytical arguments (Sec. [6]) to support this claim. Note 
that many ecological and social situations in which evolutionary dynamics take place can be 
modeled as directed or weighted networks rather than undirected and unweighted networks. 
Spreads of a type from one habitat to another may be easier than in the other direction 
because of heterogeneity in habitats and other geographical factors (see Sec. [T]for references). 
Accordingly, fixation in real contact networks may be less controlled by the values of fitness 
than in the corresponding undirected networks or well-mixed populations. 

Second, the strength of the selection pressure depends on the update rule for various net- 
works, no matter whether edges are directed or not (Sees. [Hand [5]). This conclusion extends the 
previous findings (Antal et ai, 2006; Sood et ai, 2008). For BD-D, DB-B, and DB-D, no evo- 
lutionary amplifier has been found by our exhaustive numerical analysis for small (Sec. 15. ip . 
Moreover, DB-B is strongly suppressing compared to the other update rules. In evolutionary 
game theory, DB-B (Nakamaru et al, 1998; Ohtsuki et ai, 2006; Ohtsuki and Nowak, 2006) 
and its variants (Nowak and May, 1992; Santos and Pacheco, 2005) have commonly been used. 
In games, the fitness of a type is not constant, as assumed in the present paper, but depends on 
neighbors' types. Therefore, the results obtained in this work are not immediately applicable to 
evolutionary games. However, these previous results on evolutionary games might considerably 
change under update rules less suppressing than DB-B, as has been argued (Nakamaru et ai, 
1998; Ohtsuki et al, 2006; Ohtsuki and Nowak, 2006). 

We have treated directed but unweighted networks for the sake of clarity of the analysis. 
However, we believe that our conclusions for directed networks extend to the case of weighted 
networks. This is partly supported by the analysis of the two-node weighted network presented 
in Sec. El In addition, if we quantize the edge weight and allow multiple directed edges between 
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two nodes, a weighted network can be regarded as a directed unweighted network. 

A unique point in our analysis is the use of very small networks with = 6. Small networks 
per 86 appear in many ecological contexts. For example, the number of relevant habitats may 
not be very large (Gustafson and Gardner, 1996; Tischendorf and Fahrig, 2000; Schick and 
Lindley, 2007). In addition, small directed networks are likely to be structural and functional 
building blocks of large networks (Milo et al, 2002; Itzkovitz et ai, 2003). According to the 
present analysis, the effect of the directionality of edges and that of the update rule are mostly 
consistent between small and large networks. We believe that the conclusions of the present 
paper apply to real evolutionary dynamics in populations of various scales. 
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Appendix A: Numerical procedure for small networks 

Here, we explain the methods for exactly calculating the fixation probability for strongly con- 
nected networks. Because each node takes either type A or B, there are 2^ possible states 
of the evolutionary dynamics. We define the fixation probability for a state, which may have 
multiple type-A nodes, by the probabihty that type A fixates starting from that state. The 
fixation probability for the all- A state is 1, and that for the all-B state is 0. 

Given a network, an update rule, and r, the evolutionary dynamics are equivalent to a 
nearest-neighbor random walk on the A^-dimensional hypercube comprising 2^ points. In each 
time step, the type at a node replaces the type at one of its neighbors in the original evolutionary 
dynamics. Therefore, at most one node fiips its type per unit time. In terms of the random 
walk on the hypercube, the walker moves to a neighboring state or does not move. The random 
walk continues until either the all-A state or the all-B state is reached. 
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Consider an example network shown in Fig. [71 Denote the fixation probabihty for the 
state shown in Fig. [7|by F{BAB), where the first, second, and the third arguments {A or B) 
correspond to the types at vi, V2, and V3, respectively. In a unit time, state BAB may change 
to AAB, BBB, or BAA. Otherwise, it does not change. Under BD-B, for example, the fixation 
probability satisfies 

p^BAB) = ^P(^^B) + F(BAB) ^ ^F(AAB) + FjBAA) ^ 

The first, second, and third terms in the RHS of Eq. fll4p represent the propagation of the type 
at Vi,V2, and ^3 to its neighbor, respectively. Noting the boundary conditions F{BBB) = and 
F{AAA) = 1, we can write down the other five linear equations corresponding to the single-step 
transition of F{ABB), F{BBA), F{AAB), F{ABA), and F{BAA). The fixation probabilities 
are obtained by solving the system of the six linear equations. The fixation probability starting 
from a single type-A node is given by F(r) = [F (ABB) + F (BAB) + F (BBA)] /3. Generally 
speaking, the fixation probability for an A^-node network is obtained by solving a system of 
2—2 linear equations. A standard method such as the Gauss elimination requires 
time of computation. 

For small A^, we can enumerate all the possible networks (Milo et ai, 2002; Itzkovitz et 
al, 2003). There are 1047008 strongly connected directed networks out of 1530843 weakly 
connected directed networks with A^ = 6. We cannot solve the fixation probability for A^ > 6 
due to the computational cost of enumerating the directed networks, but not due to the cost 
of solving systems of 2^ — 2 linear equations. 

Appendix B: Fixation probability for simple graphs 

In this appendix, we show detailed calculations of the fixation probability for simple graphs. 
Following a standard procedure, we map the evolutionary dynamics onto the discrete-time 
nearest-neighbor random walk on interval {0, 1, 2, ... , A^}. The position on the interval corre- 
sponds to the number of tjpe-A individuals in a network of size A^. Positions and A^ are the 
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only absorbing states. The type-A individuals increase or decrease at most by one per unit 
time. The probability that the walker at position m eventually reaches position satisfies 
the following relations: 

Fo = 0, (15) 

= amFm+l + {1 - a,n- Pm)Fm + PmFm-l, {I < m < N - 1) (16) 

Fn = 1, (17) 

for some and jSm {1 < m < N — 1). Then, the fixation probability for a single mutant is 
represented by (Moran, 1958; Ewens, 2004; Nowak, 2006): 



"T Z^m=l llm'=l a , 

Complete graph 

Consider the complete graph with nodes. The evolutionary dynamics under BD-B are 
equivalent to the Moran process, and Eq. ( |T6l) becomes 

F^ = ^ — ■ (19) 

mr + J\ — m 

Therefore, am = mr{N — m)/[{mr + N — m){N — 1)], (3m = m{N — m)/[{mr + N — m){N — 1)], 
and [3m/(y.m = 1/?", which plugged into Eq. (|T8l) leads to Eq. ([T]). DB-D and LD also yield 
(3m/otm = l/r and hence the same result as that for BD-B. 
For BD-D, we obtain 

_m {N- m)Fm+i + ^Fm N-m{N-l- m)Fm + fFm-i 

N n^^ + N-m N ^ + N - I - m ' ^ ' 

r r 

Therefore 



am r[m + r {N — 1 — m)] 
Plugging Eq. (EH) into Eq. leads to Eq. 



Pm m-l+r{N -m) ^^^^ 
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For DB-B, we obtain 

m r{m — l)Fm + (iV — m)Fm-i N — m rmFm+i + (A^ — 1 — m)F„ 
-Tm = — ; . \- 



N r{m — 1) + N — m N rm + N — \ — m 

Therefore 



j3m rm + N — 1 — m 



am r [r{m — 1) + N — m] 
Plugging Eq. ([231) into Eq. ^ leads to Eq. ^. 

Undirected cycle 

Consider the undirected cycle with nodes. Under BD-D, we obtain 

1 2 Fi + ^ N -3^ 
Fi = — Fo H — H Fi 

r 

r r 



Therefore, 



2 



Or 



r+l' (^^ = 1) 

i, (2<m<A^-2) 
%i. (m = iV-l) 



Plugging Eq. ([27D into Eq. ^ leads to Eq. ©. 
Under DB-B, we obtain 

2 rFs + Fi 1 - 3 ^ 

A* r + 1 A/ A/ 

p = "^+^^ + + Fm, (2<m< N-2) 

N r + 1 A^ r + 1 A^ rn, V - - ; 

F _ 1 p , 2 r^FV-^ + ^ AT - 3 



Therefore, 



Ml 



] {2<m<N -2) 

Clm I 2 



4^. (m = A^-l) 
Plugging Eq. (ED) into Eq. ^ leads to Eq. 
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star 

Consider the undirected star with nodes. To calculate the fixation probability for BD-B 
and LD simultaneously, it is convenient to introduce the edge weight. We assume that an edge 
from the hub to a leaf is endowed with weight unity and one from a leaf to the hub is endowed 
with weight a (Fig. [11(d)). The edge weight works as a multiplicative factor to the probability 
that this edge is chosen for reproduction, which is explained in Sec. El Denote by F^^a {Fm^B) 
{{] < m < N — 1) the fixation probability of type A when there are m type-A leaves and the 
hub is occupied by type A (B). 

BD-B for the unweighted star is equivalent to LD with selection on the birth, which is 
actually the ordinary LD, for the weighted star with a = N — 1. DB-D for the unweighted 
star is equivalent to LD with selection on the death, which is again the ordinary LD, for the 
weighted star with a = 1/{N — 1). Therefore, we calculate the fixation probability for the 
weighted star under LD. Interpreting LD as the LD with selection on the birth, we obtain 

rmaFm,A + (iV - 1 - m)aFm,B + rmFm,A + r{N - 1 - m)Fra+i,A 



m,B 



rma + {N — 1 — m)a + r{N — 1) 
(0 <m < A^-2) 

rmaFm,A + (A^ - 1 - m)aFm,B + mFm-i,B + {N - 1 - m)Fm,B 

rma + (A^ - 1 - m)a + (A^ - 1) 
{l<m< N -I) 



(32) 



(33) 



Eqs. (1521) and (15511 . respectively, lead to 



{r + a)Fr„,A = aF^^B + rFm+i,A, {0 < m < N - 2) (34) 
(ra + l)Fm,B = raFm,A + i^m-i.s. (1 < m < A^ - 1) (35) 



By combining Eqs. fl54l) and fl55l) . we obtain 

(r2a + 2r + a)F„,^-r(ra + l)F^+i,A-(r + a)F„_i,^ = 0. (1 < m < A^ - 2) (36) 
The solution to Eq. (l36l) that satisfies (r + a)Fo,A = which comes from Eq. (l34l) with 
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m = and Fq ^ = 0, and -FV-i ^ = 1 is 



r(r+a) 

ra+1 \ r{ra+l) 
r(r+a) 



ra+1 \ r(ra+l) 



Plugging Eq. ([37]) into Eq. leads to 



r(r+a) 
ra+1 



1 - 



r(ra+l) 



r(r+a) 



N-1 



ra+1 I r(ra+l) 

The fixation probability for a single mutant is represented by 



Fir) 



N-l r+a 1 ra+1 

N r(ra+l) N r(r+a) 



N 



1 _ ( ■I+^L'] 

^ ' \ra+lj 



N-2 



For BD-D on the unweighted star, we obtain 



N 



{{]<m<N -2) 



t+AT 



'-l-m -^m-l,B + in_^j^_^rn -^m,B + fT^Fm^A + (N — 1 — m)Fm,B 



N 



(37) 



(38) 



(39) 



(40) 



(41) 



(1 < m < - 1) 

Equations ( HOl) and ( HTi) . respectively, lead to 

[r{N -m)+m]Fm,A = [m + r (A^ - 1 - m)] + rF^+i,A, (42) 

(0 < m < A^-2) 

[m + l + r{N -l-m)]F^^B = [m + r (iV - 1 - m)] + i^m-i.B- (43) 

fl< m < iV- 1) 



By combining Eqs. fl42l) and fl4Hl) . we obtain 

[m + r (A^ — 1 — m)] [m — 1 + r (A^ + 1 — m)] 



m+l,A~-rm,A 



r [m + 1 + r (A^ — 1 — m)] [m — 1 + r (A^ — m) 



[l<m< N-2) 
(44) 
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which yields 



Fn-i,a — Fq^a 

^"^ [m' + r{N-l- m')] [m' - 1 + r {N + 1 - m')] 



En 

iV-2 

E 



r [m' + 1 + r (A^ — 1 — m')] [m' — 1 + r (A^ — m ^ 

N [m + r {N - 1 - m)][l + r {N - 1)] 
r"^ (A^ - 1) + 1 + r (A^ - 1 - m)] [m + r (iV - m)] 



l,A 



^0„ 



rN[l + r{N - 1)] 
(AT - l)(r- 1) 

[l + r(Ar-l)] 
(AT - l)(r - 1) 



N-2 

E 

m=0 



m + r (A^ — m) m + 1 + r (A^ — 1 — m) 



0,AJ 



1 - 



N -1 + r 



(45) 



By substituting Fi a = A^-^o,^, which is derived by setting m = in Eq. ( l42l) . and F/v_i ^ = 1 
into Eq. (H5l) . we obtain 

-^o,A = ■ 



1 



. I l+r(Af-l) . _ r-N+^N 
^ r-1 l^-^ (AT-l+r) 

By setting m = 1 in Eq. pHj) and using Fq^b = 0, we derive 

_ [i+r(Ar-2)]Ar ^ 
^^'^ - 2 + r(Ar-2) 

Therefore, the fixation probabihty F(r) = [(A^ — + Fo,a]/A^ is given by Eq. ([9]) 

For DB-B on the unweighted star, we obtain 

mF„,A + (AT - 1 - m)F„+i,A + -^^,A+(A.-i-.n)F.,, 



(46) 



(47) 



m,B 



rm+N—l—m 



N 



{0<m< N -2) 

mF^^,,B + iN-l- m)Fr.,B + '''^''-^^r:^^^-^''^" 



N 



(4J 



(49) 



(1 < m < A^- 1) 
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Eqs. fHHl) and fHUl) . respectively, lead to 



{rm + N -m)Fm,A = Em,B + {rm + N - 1 - m)F^+i,A, (50) 

(0 < m < A^- 2) 

[r{m + l) + N -l-m]Fm,B = rF^,A + (rm + - 1 - m)F„_i,B. (51) 

(1 < m < A^- 1) 



By combining Eqs. (]50l) and ( 15T|) . we obtain 

Fm+i,A - F^^A = 1 (^-,A - Fm-i,A) , (1 < m < AT - 2) (52) 

[r [m + \) + IS — \ — m\ 

which yields 

F F - V n + ^ + ^ 

^^-M-^o,A - 2^ 11 + 1) + ^ _ 1 _ (^1.^ - ^o,a) 

m=Om'=l L V ' / ' J 

(iV + ^_l)(iV_l) 
= rN-r + 1 (^^'^ - • 

By substituting Fi^^ = NFq^a/{N — 1), which is derived by setting m = in Eq. (13U|) . and 
F/v-i,yi = 1 into Eq. we obtain 

1 



1 _l_ Af+r-1 • 
rAf-r+1 



By setting m = 1 in Eq. (J5T1) and using Fq^b = 0, we derive 



(54) 



rN 

^^'^ = (Ar-l)(Ar + 2r-2)^°'^- ^^^^ 
Therefore, the fixation probability F{r) = [{N — 1)Fi,b + Fq^a\/N is given by Eq. fllOj) . 
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Figure 2: (Color online) Fixation probability for (a, b) the complete graph, (c, d) the undirected 
cycle, and (e, f) the star, under different update rules. We set N = 10 and vary r in (a, c, e). 
We set r = 4 and vary in (b, d, f). Thick solid black lines, medium solid red lines, and thin 
solid magenta lines correspond to the standard Moran process, BD-D, and DB-B, respectively. 
In (e, f), thick dashed green lines and thin dashed gray lines correspond to BD-B and DB-D, 
respectively. 
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Figure 3: Fixation probability F{A) for the strongly connected networks (gray dots) and for 
the undirected networks (black circles) of size = 6. The line F{4) = 0.7502 corresponds to 
the Moran process, (a) BD-B versus stdiT^"-), (b) BD-D versus (k), (c) DB-B versus {k), (d) 
DB-D versus stdiT"""^), and (e) LD versus (A;*"/A;°"*). 



33 




Figure 4: The network with N — 6 that yields the largest i^(4) under LD. 
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Figure 5: Fixation probability -F(4) for large undirected networks (solid lines) and directed 
networks (dashed lines) of different size N. The update rules are (a, b) BD-B, (c, d) BD-D, 
(e, f) DB-B, (g, h) DB-D, and (i, j) LD. We use the ER random graphs in (a, c, e, g, i) and 
the SF networks with p{k) oc in (b, d, f, h, j). We set (k) = 10. The Moran process yields 
F(4) ^ 0.75. 
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Figure 6: Fixation probability F{A) under DB-B for the undirected (solid line) and directed 
(dashed line) ER network with different mean degrees. We set N = 500. The Moran process 
yields F(4) ^ 0.75. 
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Figure 7: A directed network with N = 3. 
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Table 1: Correlation coefficient between F(4) and the order parameters of the network. 



Parameters 


BD-B 


BD-D 


DB-B 


DB-D 


LD 


(k) 


0.2196 


0.7719 


0.7905 


0.3791 


0.2179 


{k) (1/A;^") 


-0.3183 


-0.2639 


-0.4478 


-0.3816 


-0.1892 


(k) (1/A;°"*) 


-0.3486 


-0.6495 


-0.5415 


-0.7237 


-0.6993 


{kf/({k^n") 


0.2983 


0.3526 


0.5493 


0.4311 


0.2046 


{kf / {[k'^^'f) 


0.3392 


0.7205 


0.6395 


0.7496 


0.6921 


{kf 1 (A;™A;°"*) 


-0.4461 


-0.2715 


-0.3130 


-0.1665 


-0.6770 




-0.4640 


-0.3111 


-0.4824 


-0.3716 


-0.5048 




-0.5261 


-0.6143 


-0.5369 


-0.5978 


-0.8557 




-0.3628 


-0.4664 


-0.6233 


-0.7581 


-0.3613 


std{T"") 


-0.3790 


-0.5965 


-0.6510 


-0.8505 


-0.1935 
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Table 2: Statistics of F(4) for the networks with N = 6. 



update 


all networks 


undirected only 


rule 


(avc ± std) 


(ave ± std) 


BD-B 


0.7422 ± 0.0107 


0.7575 ± 0.0090 


BD-D 


0.5864 ± 0.0552 


0.6489 ± 0.0336 


DB-B 


0.4806 ± 0.0485 


0.5367 ± 0.0485 


DB-D 


0.6986 ± 0.0317 


0.7033 ± 0.0411 


LD 


0.7088 ± 0.0308 


0.7502 ± 
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